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Recent high-resolution A'-body CDM simulations indicate that nonsingular three-parameter models such as the Einasto profile perform 
better than the singular two-parameter models, e.g. the Navarro, Frenk and White, in fitting a wide range of dark matter haloes. While 
many of the basic properties of the Einasto profile have been discussed in previous studies, a number of analytical properties are 
still not investigated. In particular, a general analytical formula for the surface density, an important quantity that defines the lensing 
properties of a dark matter halo, is still lacking to date. To this aim, we used a Mellin integral transform formalism to derive a closed 
expression for the Einasto surface density and related properties in terms of the Fox H and Meijer G functions, which can be written 
as series expansions. This enables arbitrary-precision calculations of the surface density and the lensing properties of realistic dark 
matter halo models. Furthermore, we compared the Sersic and Einasto surface mass densities and found differences between them, 
which implies that the lensing properties for both profiles differ. 

Key words, methods: analytical, gravitational lensing: strong, gravitational lensing: weak, galaxies: clusters: general, galaxies: halos, 
cosmology: dark matter. 



1. Introduction 

The A cold dark matter (CDM) model, with (0„„ Qa) = 
(0.3, 0.7), has become the standard theory of cosmological 
structure formation. ACDM seems to agree with the ob- 
servations on cluster-sized scales (Primack 2003); however, 
on galaxy/sub-galaxy scales there appears to be a discrep- 
ancy between observations and numerical simulations. High- 
resolution observations of rotation curves, in particular of 
low surface brightness (LSB) and dark matter dominated 
dwarf galaxies (deBloketal. 2001; van den Bosch & S waters 
2001; Swatersetal. 2003; Weldrake et al. 2003; Donatoetal. 
2004; Gentile et al. 2005; Simon et al. 2005; Gentile et al. 2007; 
Banerjee et al. 2010) favour density profiles with a flat central 
core (e.g. Burkert 1995; Salucci & Burkert 2000; Gentile et al. 
2004; Li & Chen 2009). In contrast, A^-body (dark matter only) 
CDM simulations predict galactic density profiles that are too 
high in the centre (e.g. Navarro, Frenk, & White (1996, 1997) 
(NFW); Moore et al. 1999). This discrepancy is called the cusp- 
core problem; a complete review can be found in de Blok (2010). 

Gravitational lensing is one of the most powerful tools in 
observational cosmology for probing the matter distribution 
of galaxies and clusters in the strong regime (Kochanek et al. 
1989; Wambsganss & Paczynski 1994; Bartelmann 1996; 
Chae et al. 1998; Kochanek et al. 2000; Keeton & Madau 2001; 
Sandetal. 2002; Keeton 2002, 2003; Keeton & Zabludoff 
2004; Limousin et al. 2008; Anguita et al. 2009; Zitrin et al. 
2011a,b) and the weak regime (Kaiser & Squires 1993; 
Mellier 1999; Bartelmann & Schneider 2001; Hoekstra et al. 
2004; Clowe et al. 2006; Mahdavi et al. 2007; Jee et al. 2009; 
Huang et al. 201 1). Comparing these observations to theoretical 
models provides key information to help resolve the cusp-core 
problem. Evidently, one must use the most accurate density pro- 



file to obtain the best fit to observational data from strong- and 
weak-lensing studies. 

Recently, A^-body CDM simulations (Navarro et al. 2004; 
Merrittetal. 2006; Gaoetal. 2008; Hayashi & White 2008; 
Stadel et al. 2009; Navarro et al. 2010) have found that cer- 
tain three-parameter profiles provide an excellent fit to a wide 
range of dark matter haloes. One of these is the Einasto (1965) 
profile, a three-dimensional version of the two-dimensional 
Sersic (1968) profile used to describe the surface brightness 
of early-type galaxies and the bulges of spiral galaxies (e.g. 
Daviesetal. 1988; Caon et al. 1993; D'Onofrio et al. 1994; 
Cellone et al. 1994; Andredakis et al. 1995; Prugniel & Simien 
1997; Mollenhoff'&Heidt 2001; Graham & Guzman 2003; 
Graham et al. 2006; Gadotti 2009). The Sersic profile can be 
written as: 



l/m 



1 



(1) 



where R is the distance in the sky plane, m the Sersic index, 
T is the mass-to-light ratio, /g is the luminosity density at the 
effective radius R^., and b„, is a dimensionless function of m 
that can be determined from the condition that the luminosity 
inside R^. equals half of the total luminosity. Numerical solu- 
tions for b„, are given by Ciotti (1991), Moriondo et al. (1998), 
Prugniel & Simien (1997) and an asymptotic expansion b„, = 
2m - 1/3 -I- 4/405m + O (m-) using analytical methods was ob- 
tained by Ciotti & Bertin (1999). 

The Einasto profile (model) is characterised by a power-law 
logarithmic slope, 

dlnp 



y(r) = - 



dlnr 



(r) oc r'/". 



(2) 



with «, which we call the Einasto index, a positive number defin- 
ing the steepness of the power law. Integrating leads to the gen- 
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eral density profile 
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(3) 



where represents the radius of the sphere that contains half 
of the total mass, ps is the mass density at r = r^, and d„ is a 
numerical constant that ensures that rs is indeed the half-mass 
radius. In the context of dark matter haloes, the density can also 
be expressed as 



P W = P-2 exp \ -In 



i/« 

— I -1 



(4) 



where p_2 and r_2 are the density and radius at which p (r) cx r"^ . 
In the remainder of this paper, we will use yet another, equivalent 
version, 

p(r) = po exp 



-(3 



(5) 



where we introduced the central density 



Po = Ps e''" = P-2 e^'\ 



and the scale length 



h - — - 



(6) 



(7) 



If a model is to describe real galactic systems, several condi- 
tions must be imposed on the model description functions, as 
discussed by Einasto (1969a). When a model is constructed, an 
initial descriptive function is chosen; the most practical choice is 
the density profile, because the main descriptive functions (cu- 
mulative mass profile, gravitational potential, surface mass den- 
sity) are integrals of the density profile. Furthermore, a physi- 
cal model has to satisfy several conditions: i) the density pro- 
file must be non-negative and finite; ii) the density must be a 
smoothly decreasing function that approaches zero at large radii; 
iii) some moments of the mass function must be finite, in partic- 
ular moments that define the central gravitational potential, the 
total mass, and the effective radius of the system; and iv) the de- 
scriptive functions must not have jump discontinuities. Einasto 
(1969a) presented several families of valid descriptive functions, 
among which the so-called Einasto profile is a special case that 
agrees best with observations. 

The Einasto profile was used by Einasto (1969b) to obtain a 
model of M31. Later, this model was applied by Einasto (1974) 
to several nearby galaxies, including M32, M87, Fornax and 
Sculptor dwarfs, and also M3 1 and the Milky Way. These mod- 
els were multi-component ones, and each component has its own 
parameter set {po, h, n } representing certain physically homo- 
geneous stellar populations. 

Navarro et al. (2004) found that for haloes with masses 
from dwarfs to clusters, 4.54 < « < 8.33 with an average 
value of n = 5.88. Hayashi & White (2008) and Gaoetal. 
(2008) found that n tends to decrease with mass and redshift, 
with n ~ 5.88 for galaxy-sized and n ~ 4.35 for cluster- 
sized haloes in the Millennium Run (MR) (Springel et al. 2005). 
Navarro etal. (2010) obtained similar results for galaxy-sized 
haloes in the Aquarius simulation (Springel et al. 2008). Also, 
Gaoetal. (2008) showed that n ~ 3.33 for the most massive 
haloes of MS. Chemin et al. (201 1) modelled the rotation curves 
of a spiral galaxies subsample from THINGS (The HI Nearby 
Galaxy Survey, (de Blok et al. 2008)), using the Einasto profile. 



and found that n tends to have lower values by a factor of two 
or more compared to those predicted by A^-body simulations. 
Dhar & Williams (201 1) fitted the surface brightness profiles of 
a large sample of elliptical galaxies in the Virgo cluster, using a 
multi-component Einasto profile consisting of two or three com- 
ponents for each galaxy with different n for each individual com- 
ponent. For the central component, they found values of n < 1 
for the most massive and shallow-cusp galaxies and « < 2 for 
steep-cusp galaxies, while 5 < « < 8 for the outer component of 
massive galaxies. These values for the outer component are con- 
sistent with the results from A^-body simulations for Einasto dark 
matter haloes, and the authors argued that Einasto components 
with an n in this range could be dark matter dominated. 

In the light of its increasing popularity to describe the den- 
sity of simulated dark matter haloes, a detailed investigation 
of the properties of the Einasto model is of paramount impor- 
tance. Some aspects of the Einasto model have been presented 
by several authors (Mamon & Lokas 2005; Cardone et al. 2005; 
Merritt et al. 2006; Dhar & Williams 2010). The most complete 
study of the properties of the Einasto model is the work by 
Cardone et al. (2005), who provided a set of analytical expres- 
sions for quantities such as the mass profile and gravitational po- 
tential and discussed the dynamical structure for both isotropic 
and anisotropic cases. Nevertheless, the Einasto model has not 
been studied analytically as extensively as the Sersic models, and 
several properties still have to be further investigated in more de- 
tail. One area where progress is still to be made is on the value 
of the dimensionless parameter d,,- 

The most important lacuna, however, concerns the surface 
density on the plane of the sky, an important quantity that de- 
fines the lensing properties of a dark matter halo. The sur- 
face mass density has an importance in theoretical predictions 
and observations. On the observational side, based mostly on 
the mass decomposition of rotation curves using cored haloes, 
various studies (Kormendy & Freeman 2004; Spano et al. 2008; 
Donato et al. 2009; Gentile et al. 2009) have shown that the 
product of the central density po and the core radius ro is con- 
sistent with a universal value, independent of galaxy mass. The 
product poro is directly proportional to the average surface den- 
sity within ro, and to the gravitational acceleration of dark mat- 
ter at ro. In view of a detailed comparison of these results with 
the outcome of the most recent dark matter simulations, it is 
therefore important to study the analytical properties of the sur- 
face density distribution of Einasto haloes. Recent works study- 
ing theoretical and observational aspects of dark matter surface 
densities include Boyarsky et al. (2010), Walker et al. (2010), 
Napolitano et al. (2010), and Cardone etal. (2011): efforts are 
being made to confirm or call into question the universality of 
the dark matter surface density within one dark halo scale length. 

Cardone et al. (2005) showed that the surface density of the 
Einasto model cannot be expressed in terms of elementary func- 
tions and discussed the general properties using numerical inte- 
gration. Dhar & Williams (2010) presented an analytical approx- 
imation for the Einasto surface brightness profile and demon- 
strated that this surface density profile is not Sersic-like. A 
general analytical formula, which would enable an arbitrary- 
precision calculation and an analytical study of the asymptotic 
behaviour, is still lacking to date. 

The most recent studies in strong and weak lens- 
ing (e.g. Donnarumma et al. 2011; Sereno & Umetsu 2011; 
Morandi et al. 201 1; Umetsu et al. 201 1; Oguri et al. 201 1) use 
the NEW profile or its generalization (Zhao 1996; Jing & Suto 
2000) to model the dark matter halo instead of the Einasto pro- 
file. One of the factors that limits its adoption is the absence of 
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analytical formulas for its lensing properties. Therefore, its ap- 
plication in lensing studies is not as wide as for other profiles. 
A complete general set of analytical formulas for the lensing 
properties of the Einasto profile would help to increase its use in 
lensing studies. 

In this paper, we extend the analytical study of the Einasto 
model using some of the techniques also employed in the 
extensive literature on the analytical properties of the Sersic 
model (e.g. Ciotti 1991; Ciotti & Lanzoni 1997; Ciotti & Bertin 
1999; Trujillo et al. 2001; Mazure & Capelato 2002; Cardone 
2004; Graham & Driver 2005; EKasddttir & Moller 2007; 
Baes & Gentile 201 1; Baes & van Hese 201 1). In Section 2, we 
discuss an analytical expansion for the dimensionless parame- 
ter d„. In Section 3, we derive an analytical expression for the 
Einasto surface mass density in terms of the Fox H function, 
using the Mellin transform-method. We then use the result for 
the projected surface mass density to calculate the cumulative 
surface mass, lens equation, deflection angle, deflection poten- 
tial, magnification, shear and the critical curves for a spherically 
symmetric lens described by the Einasto profile in terms of this 
function. We also calculate explicit series expansions for the sur- 
face mass density and all the lensing properties. In Section 4, 
we derive some special cases for the lensing properties for inte- 
ger and half-integer values of n in terms of the Meijer G func- 
tion. Furthermore, we compare the Einasto and Sersic surface 
mass densities for the same values of their respective indices. 
Finally, in Section 5 we discuss the asymptotic behaviour of the 
surface mass density and cumulative surface mass at small and 
large radii. We present our conclusions in Section 6. 

2. Spatial properties 

The total mass of an Einasto model with a mass density profile 
given by (5), is 

M = 47rpo/^^n^(3n). (8) 

If we use this formula to replace the central density po by the to- 
tal mass M as a parameter in the definition of the Einasto models, 
we obtain 

p(n = ; e 

^ 47r/i3«r(3«) 

where we have introduced the reduced radius 

(d„r r 
s = . 

rs 

At small radii, the density profile behaves as 
M 



(9) 



(10) 



p{r) 



47ThhiT{3n) 



(l-.'/" + ---). 



(11) 



The cumulative mass profile M(r) of a spherical mass distribu- 
tion is found through the equation 



M{r) = An 
For the Einasto model, we find 



r p{r')r'^dr'. 
Jo 



M{r) = M 



1 



r(3«, s^'") 



r(3«) 



where T(a, x) is the incomplete gamma function. 



T(a, x) 



(12) 



(13) 



(14) 



Given that this is the radius of the sphere that encloses half of 
the total mass, we find that d„ is the solution of the equation 



2r(3n,c/„) = r(3«). 



(15) 



This equation cannot be solved in a closed form; one option is to 
solve it numerically or use interpolation formulae. For example, 
Merritt et al. (2006), quoting G. A. Mamon (private communi- 
cation), propose d„ ^ 3n - ^ + 0.0079/«, but do not mention 
the origin of this approximation or quote its accuracy. An el- 
egant way to determine an approximation for d„ builds on the 
work by Ciotti & Bertin (1999), who used an asymptotic expan- 
sion method to solve the general relation r{a, x) - r{a)/2. If we 
apply the resulting expansion to our problem, we obtain for the 
parameter d„ from the Einasto model 



1 8 
d„ ~ 3n - — + . , 
3 1215 « 

1048 

31000725 «3 ■ 



184 



229635 «2 

17557576 ^ . 
1242974068875 \ « ' 



(16) 



The gravitational potential of a spherically symmetric mass dis- 
tribution p{r) can be found through the formula 



>I'(r) = 47tG 
or equivalently 



p(r')r'^dr' + J pir')r'dr' 



r 



For the Einasto model, we find 
GM 



M(r') dr' 



r(3n) r(3«) 

The Einasto model has a finite potential well, given by 

^ _ GM T{2n) 
° h T(3n)' 

At small radii, the potential decreases parabolically 
T(2«) 



h 



r(3«) 6nr(3n) 



(17) 



(18) 



(19) 



(20) 



(21) 



which is not surprising, given the finite density core of the 
Einasto density profile. At large radii, we obtain a Keplerian fall- 
off, 

^(0- — + ■••, (22) 
r 

which is also expected for a model with a finite total mass. 



3. Lensing properties 

3.1. Surface mass density and cumulative surface mass 
density 

The surface mass density of a spherically symmetric lens is given 
by integrating along the line of sight of the 3D density profile 



2(^) 



-f 

•J —a 



P(€,r) dz. 



(23) 
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where ^ is the radius measured from the centre of the lens and 
r — ^J^~ + z~. This expression can also be written as an Abel 
transform (Binney & Tremaine 1987) 



I 



p (r) r dr 



Inserting eq. (5) into the above expression, we obtain 



S (x) = 2 po 



r 



e sas 



(24) 



(25) 



where we have introduced the quantities x = ^/h and s - r/h. 
As discussed by Cardone et al. (2005) and Dhar & Williams 

(2010) , the integral (25) cannot be expressed in terms of ele- 
mentary or even the most regular functions for all the values of 
n. Only the central surface mass density can be evaluated analyt- 
ically as 

2:(0) = 2«po/^^(«). (26) 

This situation is very similar to the deprojection of the Sersic 
surface brightness profile (1). Indeed, the deprojection formula 
leads to an integral that is an inverse Abel transform, and also in 
this case it was long thought that no further analytical progress 
was possible. However, Mazure & Capelato (2002) used the 
computer package Matheraatica to demonstrate that it is possi- 
ble to write the Sersic luminosity density in terms of the Meijer 
G function for all integer Sersic indices m. Baes & Gentile 

(201 1) and Baes & van Hese (201 1) took this analysis one step 
further and showed that the deprojection of the Sersic surface 
brightness profile for general values of m can be solved elegantly 
using Mellin integral transforms and gives rise to a Mellin- 
Barnes integral. The result is that the Sersic luminosity density 
can be written compactly in terms of a Fox H function, which 
reduces to a Meijer G function for all rational values of m. 

The obvious similarity between these two cases invites us 
to apply the same Mellin transform technique (Marichev 1982; 
Adamchick 1996; Fikioris 2007) to the integral (25). The basic 
idea behind this technique is that any definite integral 



f(z) 



f 

Jo 



g(t,z)dt. 



can be written as 



(27) 



(28) 



This expression is exactly a Mellin convolution of two functions 
/i and /2. Now we can apply the Mellin convolution theorem, 
which states that the Mellin transform of a Mellin convolution 
is equal to the products of the Mellin transforms of the origi- 
nal functions. As a result, the definite integral (27) can be writ- 
ten as the inverse Mellin transform of the product of the Mellin 
transforms of fi and /2 . The Mellin transform and its inverse are 
defined as 



f(z)z"-'dz, 

a 

^-/(z)^f(z)^^.j^<P(u)z-"du, 



(29) 
(30) 



with X. a line integral over a vertical line in the complex plane. 
This implies that 



fiz) 



= -^-^ f m /, (m) to (u) z'" du. (31) 
27ri 



It now turns out that the integral (31) is a Mellin-Barnes integral 
for large classes of the functions f\ and /2, and that this integral 
can be evaluated as a Fox H function or a Meijer G function in 
many cases. 

We can immediately apply this formalism to the integral 
(25), with z = 1 and 



/i(f) = 2po/je-'""f\ 



(32) 



and 



, if < f < x^i, 

m)^\ Vl - X2f2 (33) 

[ else. 
The Mellin transforms of these functions are readily calculated 
%,(m) = 2po/!nr[n(2-HM)], (34) 



Combining these results, we obtain 
2(x) = 2n ^p^hx 

1 ^ r(2n3')r(-i+3;) 



(35) 



r(3') 



[x^p'dy. (36) 



This Mellin-Barnes integral may be recognized as a Fox H func- 
tion, which is generally defined as the inverse Mellin transform 
of a product of gamma-functions. 



Tjmji 



(a. A) 




(b,B) 


z 



i_ r 



n™ 1 nbj + Bjs) Y\']_ J(l-aj-Ajs) 
ni-bj-Bjs)Y\']__,^^J(,aj+Ajs) 



z'-'ds. 



(37) 



Using this definition, we can write the surface density of the 
Einasto model in the following compact form. 



2 (x) = 2n -s/tt po h x H 



(0, 1) 
(0, 2«),(-i, 1) 



x2 



(38) 



The Fox H function is a general analytical function that is 
becoming more and more used both in mathematics and ap- 
plied sciences; in fact, it is scheduled for inclusion in the 
Mathematica numerical library. While not the most mainstream 
special function, the Fox H is an extremely flexible function 
that contains very broad classes of elementary and special func- 
tions as particular cases. It is indeed a very powerful tool for 
analytical work. It has many general properties that allow one 
to manipulate expressions to equivalent forms, reduce the or- 
der for certain values of the parameters, etc. For details on its 
many useful properties, we refer the interested reader to Mathai 
(1978), Kilbas & Saigo (2004), or Mathai et al. (2009) and the 
references therein. 

As an illustration of the power of the Fox H function as an 
analytical tool, and as a sanity check on the formula, we can 
calculate the total mass of the Einasto model by integrating the 
surface mass density over the plane of the sky. 



2;(^)^d^ = 2«7r^/-po/i 
Jo 



(0, 2«),(-i 1) 



t^'^dt. (39) 
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obtain after some algebra 



This integral can be calculated by setting ^ = 5 and a = 1 in 
equation (2.8) in Mathai et al. (2009). We immediately find 



M^lnpo n ^''^"•'3^''^^ = AnpQ Ir' n Y (3n) , (40) 



in agreement with equation (8). 

An important quantity for gravitational lensing studies is the 
cumulative surface mass density, i.e. the total mass contained in 
a infinite cylinder with radius ^, 



M{^) = In 



(41) 



We find 



M(x) = Inn^'^poh 



3 3 
X 



X H. 



i-j, 1),(0, 1) 
(0, 2n),(-l !),(-!, 1) 



(42) 



As another demonstration of the usefulness of the Fox H func- 
tion, we can apply the residue theorem to the contour inte- 
gral (37), and obtain explicit series expansions. The general 
procedure can be found in Kilbas & Saigo (1999), and a spe- 
cific application to the deprojected Sersic model is given in 
Baes & van Hese (201 1). The analysis for the projected Einasto 
profile is completely analogous: again, the form of the series ex- 
pansion depends on the multiplicity of the poles of the gamma 
functions Fibj + B js). For both i;(x) and M{x), the poles of these 
gamma functions are -kijln and 1/2 - k2, with ki and ki any 
natural number. We encounter two cases: 

Case 1: if « is either non-rational or a rational number pjq 
with an even denominator (and p, q coprime), all poles are sim- 
ple and the expansion is a power series. 



S(jic) — In VtFpo h 



£ ^(-^-1;) (_i)^ ^^/« 



k\ In 
Y V(n-lnk) {-Xf 



,1k 



(43) 



and 



M(x) = InTv'l-pah^ 



r(-i-|;) (-1)^ x"!"^^ 



-^Tin-Ink) {-If 



k- 



2n 



2k+2 



(44) 



2](x) = 2n po h 



r(«) nn — z.r 

1^ ^ h rfi-- 



k:i A-i) 



2n 



Y{ n-2nk) {-If 
k. ^{h-k) ^' ' 

(k+ko) mod qi=0 



. , {-l)P{2k)\ lx\^k r 

-^P^'I. {2nk-nV.klkl (2) ^H^y 



(k+ko) mod q=0 

1 



+ilf{k+l) + n iff{2nk - n) - iff{2k - 1) 

2k 



(45) 



and 



M{x) = 2nn^'^poh^ x^ 



y {-Ifx"/"-^ 

^ r(-±) 



k=l 
k mod pi^{) 



k\ 2n 



T{n) 



n) sr-i i(n - 

^ ^ "r(i-fc) ik+\f 

ik+ko) mod qi^iy" 



r{n-2nk) {-If 2k 



+ 2npoh^x^ ^ 



{-l)"{2k) 



, „ {2nk-nfk\{k+ 1) 

k=0 
(k+ko) mod q=0 



1 



1 



+ i//{k + 2) + n i//{2nk - n) - i//{2k - 1) 



(46) 



2k 2k + 2 
with iff{k) the digamma function. 

3.2. Lens equation and deflection angle 

In the thin lens approximation, the lens equation for axially sym- 
metric lens is 

?1^^^-Dls&, (47) 

where the quantities rj and ^ are the physical positions of the 
source in the source plane and an image in the image plane, re- 
spectively, a is the deflection angle, and Dl, Ds and Dig are the 
angular distances from observer to lens, from observer to source, 
and from lens to source, respectively. 

With the dimensionless positions y = Dirj/Ds h and x - 
^/li, and dimensionless a - Disa/Ds ^, the lens equation 
reduces to 

y - X- a{x). (48) 

The deflection angle for a spherical symmetric lens is 
(Schneider et al. 1992) 

a{x)^- x'-^^djc' = - xK{x)dx\ (49) 

X Jo ^crit X Jo 



where 



S(x) 



(50) 



Case 2: if « is integer or a rational number p/q with an odd 
denominator, some poles are of second order, and the expansion 
is a logarithmic-power series. If we define ^0 = then we 



is the convergence and Eciit is the critical surface mass density 
defined by 

2^crit = -j -prp: p; > (51) 



5 



E. Retana-Montenegro et al.: Analytical properties of Einasto dark matter haloes 



where c is the speed of Ught, and G is the gravitational constant. 
Evidently, the deflection angle is related to the integrated mass 
as 

M{x) 



a{x) = 



n h- Sni 



(52) 



Introducing the central convergence, /Cc, a parameter that deter- 
mines the lensing properties of the Einasto profile. 



Z(0) 2pQhnT(n) 



(53) 



we can write a (x) in the form 



a{x) 



r(n) 



V^r 2 „2,i 



(-i 1),(0, 1) 
(0, 2n),(-{, 1),(-|, 



1) 



(54) 



with a completely analogous series expansion as for M(x). For 
a spherically symmetric lens that is capable of forming multi- 
ple images of the source, one sufficient condition is > 1 
(Schneider et al. 1992). In the case < 1, only one image of 
the source is formed. In addition to the condition /Cc > 1, multi- 
ple images are produced only if | y |< y^t (Li & Ostriker 2002), 
where jciit is the maximum value of y when x < or its mini- 
mum when X > 0. For singular profiles such as an NEW profile, 
the central convergence always is divergent, hence the condition 
/Cc > 1 is always met, this implies that an NEW profile is capable 
of forming multiple images for any mass. Nonsingular profiles, 
such as the Einasto profile, are not capable of forming multiple 
images for any mass. Instead, the condition > I sets a mini- 
mum value for lens mass required to form multiple images. 

3.3. Deflection potential 

The deflection potential i/r (x) for a spherically symmetric lens is 
given by (Schneider et al. 1992) 



i^(x) 



2^ X k{x') In^-^j dx'. (55) 



Inserting eq. (38) into (55), we obtain again a result that can be 
re-expressed as a Eox H function 



i^(x) 



2r(«) 



X H 



2.2 



(-i, l),(-i, 1),(0, 1) 
(0, 2«),(-i, 1),(-|, 1),(-|, 1) 



x2 



Eor i^(x) we have the same cases of expansions, a power series 
for simple poles in case 1, 



iA(x) 



2r(«) 



S^j (-1)*^ X' 



.k/rn-3 



k\ 3« + k 



Y{n-2nk) x 



S r(i-^) (^+1)!^ + 1 



(57) 



and a logarithmic -power series for multiple poles in case 2, 



'/'(x) — JS^. x^ 



-z 



rH-4) (-1)' ^ 



k/n+l 



k\ 3n+k 



2r(«) ^ r{-i) 

r(n) -s^ r(n-2nk) (-if x^* 



Vi {ri rU-k) (k+iy.k + i 

(k+ko) mod qtO 

Kc 2 y (-iy(2ky. /xv 
2r(«)^ f-^ {2nk-ny.(k+iy.(k+iy.\2) 

(k+ko) mod q=0 

- — + >ff(k + '2) + n ilf(2nk - n) - il/(2k - 1) 
\2J 2k 



(58) 



3.4. Magnification, shear and the critical curves 

Gravitational lensing effect preserves the surface brightness 
but it causes variations in shape and solid angle of the 
source. Therefore, the source luminosity is amplified by 
(Schneider et al. 1992) 

^ = n 72 1' (59) 

(1 - kY - 

where k is the convergence and j - y{x) is the shear. The am- 
plification fi has two contributions: one from the convergence, 
which describes an isotropic focusing of light rays in the lens 
plane, and the other which describes an anisotropic focusing 
caused by the tidal gravitational forces acting on the light rays, 
described by the shear Eor a spherical symmetric lens, the shear 
is given by (Miralda-Escude 1991) 



y(x) = 



where the average surface mass density within x is 
2 



2(x) 



Jo 



2;(x') dx'. 



(60) 



(61) 



The magnification of the Einasto profile can be found combining 
eqs. (38), (59), (60) and (61). In the calculation of £ (x), we use 
again the Mellin-Barnes integral representation (36), and inte- 
grate it to obtain an expression in terms of Eox H function. Thus 
we get 



ti = [{l--K){l+K-2K)]- 



where 



k(x) 



r(«) 



(0, 1) 
(0, 2n),(-i, 



1) 



(56) and 



k{x) = 



r(«) 2,3 



(0. 



2n),i- 



1),(0, 1) 

i IX (- 



1) 



(62) 



(63) 



(64) 



The magnification may be divergent for some image positions. 
The loci of the diverging magnification in the image plane are 
called the critical curves. We see from eq. (62) that the Einasto 
profile has one pair of critical curves. The first curve, 1-^ = 0, 
is the tangential critical curve, which corresponds to an Einstein 
Ring with a radius, called the Einstein radius. The second curve, 
1 + ^ - 2a: = 0, is the radial critical curve, which also defines 
a ring and its corresponding radius. In both cases the equations 
must be solved numerically. 
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4. Integer and half-integer values of n 



The expressions for the Einasto surface mass density (38) and its 
lensing properties (42), (54), (56), (62), (63) and (64) in terms 
of the Fox H function for case « is a rational number, can be 
reduced to a Meijer G function. 



2m n 



■z^'ds. (65) 



Meijer G function numerical routines had been implemented 
in computer algebraic systems, such as the commercial Maple, 
Mathematica and the free open-source Sage and mpmath li- 
brary; in contrast, there is no Fox H function implementation 
available. 

For case n integer or half-integer, we can write I.ix) as a 
Meijer G function. Inserting the Gauss multipUcation formula 
(Abramowitz & Stegun 1970) 

2h-1 

T(2ny) = (2n)-^2"v (^2n)i-" r(y) F + yj , (66) 

into eq. (36) and comparing with the definition (65), we obtain 
an expression for the surface mass density of the Einasto profile 
in terms of the Meijer G function 



2(x) 



yfnpoh 2«,o 
(2;r)«-i 





X" 


b 


{2nfA 



where ft is a vector of size 2n given by 

2« - 1 



! 2n 2n 



2n 



(67a) 



(67b) 



Integrating the above expression for surface mass density ac- 
cording to (41) and using the integral properties of the Meijer G 
function (eq. 07.34.21.0003.01 on the Wolfram Functions Site'), 
we obtain 



M(x) 



2(27r)"-2 





1 1 

x~ 


A 


{2nfA 



(68) 



Inserting eq. (67) into (49) and performing the integral of the 
Meijer G function, we find 



a{x) 



2(27r)''-i Vnr(«) 



2 ,^2;i,I 
<5l,2„+l 





X 




{2nf\ 



(69) 



We obtain the deflection potential in terms of the Meijer G func- 
tion combining eqs. (55), (67) and again we integrate, applying 
the Meijer G function properties. We obtain 



4(27r)"-i Vnr(n) 



3 ^2«,2 



1 1 


X^ 1 


2' 2, 
2' 2 


{2nf"\ 



(70) 



The convergence can be simply found dividing eq. (67) by 



k{x) 



-,2/7,0 



X G' 

2(27r)«-i Vnr(«)^ 





x2 1 


b 


(2«)2"J 



(71) 



the average convergence can be calculated by inserting eq. (67) 
into eq. (61), and then integrating, we find 



k{x) 



-,2(1,1 



X G^ 

2(27r)"-i yfJirin)^ '-2"+' 





x" 1 


A 


{2nfA 



(72) 



http://functions.wolfram.com/HypergeometricFunctions/MeijerG/ 



4.1. Simple cases: n - 1 andn - \ 

The Einasto profile corresponding to n = 1 is a simple exponen- 
tial model, characterised by the density profile 



p{r) =po exp|-(^j 



(73) 



the surface density profile and lensing properties are readily 
found by setting n = 1 in eqs. (67), (68), (69), (70), (71) and 
(72) 



I.{x) - pQhxG^^l 



M{x) = npoh^ x^ Gjj 



1 I 



a{x)^ —X Gj_3 



2' 2' 2 
_ !_ 

2 



I t \ '''c 3 ^2,2 
i/r (jc) = — JC G. • 



2,4 



i _i -1 

2' 2' 2 

_i _i 

2' 2 



!_ 1 

2' 



3 _3 
2 ' ~i 

4 



-I \ c „2,1 



i _i 

2' 2 

_\ 

1 _l'_3 
2' 2' 2 



(74) 
(75) 
(76) 
(77) 
(78) 
(79) 



The specific Meijer G functions in these expressions can be re- 
duced to more standard special functions. We obtain the equiva- 
lent expressions 



S (x) — 2po h X K\ {x) , 

^2 



Mix) = 87rpo/j 



1 



■K2{x) 



4a- 

a (x) = 

X 



l-^Kiix) 



ip (x) = AKc 



In 



[^A-KAx).K,ix).y-\ 



K (x) — K^xKl (x) 

4/f, 



Kix)^ — 

X^ 



x2 

l-y^:2(-^) 



(80) 
(81) 

(82) 

(83) 
(84) 
(85) 



where Ky{x) is the modified Bessel function of the second kind 
of order v, and y x 0.57721566 the Euler-Mascheroni constant. 
The expressions (80), (81), (82), (83), (84) and (85) can also be 
readily calculated by inserting the density profile (73) into the 
formula (24) and later carrying out all corresponding calcula- 
tions for the lensing properties. 

For n - I, the density profile (4) becomes a Gaussian, 



p(r) =po exp 



(86) 
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the surface density and lensing properties can be found through 
substitution of n = ^ in eqs. (67), (69), (68), (70), (71) and (72). 



The Sersic Profile 



M (x) ^ n^'^ po x^G\i 



a{x) - K^x Gj 2 



'2 

1 h 
2' 2 



I 

1 2 3 




x~ 


_ 2' 2 





_i _i 

2' 2 
i -1 -1 
"2' 2' 2 



/f (jc) = X G^;" 



1 

1 2 3 


X^ 


2' 2 





(87) 
(88) 
(89) 
(90) 
(91) 
(92) 



These Meijer G function can also be written in terms of elemen- 
tary and special functions. 



S(ji:) = yfnpQh& , 


(93) 


M(x) ^TT' Impair' (l -e"''), 


(94) 


a{x)^- (l -e-'"'), 
X ^ ' 


(95) 


iA(^) = y [\n[x'-) + 'E,{x^) + y\, 


(96) 


K (x) - Kc e" '" , 


(97) 


,(.) = |(l-e-'-^). 


(98) 



with Ey{x) the exponential integral of order v. The above results 
can be easily checked by substituting eq. (86) into recipes (24), 
(41), (49), (50), (55) and (61). 

4.2. Profile comparison 

We compared the Einasto and Sersic surface mass densities for 
the same values of the Sersic index m and Einasto index «, in- 
cluding the exponential and Gaussian cases. For this comparison 
we used the Mathematica implementation of the Meijer G func- 
tion and eq. (67). Figure 1 shows "LsiR) for different values of m, 
while T/e and Rf. are held fixed; Figure 2 displays Z(x) for differ- 
ent values of «, while por^ and are held fixed. In both figures, 
it can be clearly seen that the respective index is very important 
in determining the overall behaviour of the curves. 

The Sersic profile is characterised by a steeper central core 
and extended external wing for higher values of the Sersic index 
m. For low values of m the central core is flatter and the exter- 
nal wing is sharply truncated. The Einasto profile has a similar 
behaviour, with the difference that the external wings are more 
spread out. Also in the inner region for both profiles with low 
values of the respectively index we obtain higher values of Es 
and E. Additionally, the Einasto profile has higher values of the 
central surface mass density than the Sersic profile, comparing 
both profiles for the same index. However, the Einasto profile 
seems to be less sensitive to the value of the surface mass den- 
sity for a given n and radius in the inner region than the Sersic 
profile. It is in this region that the lensing effect is more impor- 
tant and the surface mass density characteristics determine the 
lensing properties of the respective profiles. 

Given these differences between the two profiles, we clearly 
see that the lensing properties of Sersic and Einasto profiles also 
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Fig. 1. Sersic profile, in dimensionless units, where T/e and 
are held fixed for different values of the Sersic index m. 

The Projected Einasto Profile 
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Fig. 2. Projected Einasto profile, where p^rs and are held fixed 
for different values of the Einasto index «. 



differ, which agrees with previous work of Cardone et al. (2005) 
and Dhar& Williams (2010). Studies of the lensing properties 
of the Sersic profile have been made by Cardone (2004) and 
Eliasdottir & Moller (2007). 



5. Asymptotic behaviour 

The series expansions allow us to directly investigate the be- 
haviour of Einasto models at small radii. It follows that the cen- 
tral asymptotic behaviour of the surface density 2(x) depends on 
the value of «. If « < 1, we find the following expansion at small 
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radii (x «c 1) 

~ pah[ 2T(n + 1) - r(l - «) ] . 
\in— 1, then the expansion has the form 

l.{x)~ poh[2 + [2\n[\)-\)x']. 
Finally, if n > 1 , the central surface density behaves as 



I(x) ~ pqH 



2T{n + 1) ■ 



n + 1 



(99) 



(100) 



(101) 



The behaviour of the cumulative surface mass M{x) and deflec- 
tion potential if/ (x) are more straightforward. At small radii, we 
find the asymptotic expressions 



M(x) ~ 2npQlr'Y{n + l)x% 



and 



(102) 



(103) 



The slope for M{x) is not unexpected, given that the Einasto 
models have a finite central surface mass density. The asymp- 
totic behaviour of the Fox H function at large radii is described 
in Kilbas & Saigo (1999). We obtain the following expansions 
(x» 1) 



S(x) ~ ^innpohe. x' , 
M(x) -AnpQh^nTOn) - 2(2mi)^'^ po x^- 
and 



(104) 
(105) 



iA(x) 



,r(3«) 
r(n) 



In 



(^)-«(A(3«) 



+ 1 



■<Jn (2n)^^^ _.j/„ 

+ 77—- X 2„ . 

r(n) 



(106) 



6. Summary and conclusions 



We studied the spatial and lensing properties of the Einasto pro- 
file by analytical means. For the spatial properties we applied 
the method used by Ciotti & Bertin (1999) to derive an analyti- 
cal expansion for the dimensionless parameter d„ of the Einasto 
model. We also derived analytical expressions for the cumula- 
tive mass profile M(r) and the gravitational potential *F(r). For 
the lensing properties, we used the Mellin integral transform 
formalism to derive closed, analytical expressions for the sur- 
face mass density E (x), cumulative surface mass M (x), deflec- 
tion angle a (x), deflection potential 1^ (x), magnification ju (x) 
and shear 7 (x) . For general values of the Einasto index «, these 
are expressed in terms of the Fox H function. Using the prop- 
erties of the Fox H function, we calculated explicit power and 
logarithmic-power expansions for these lensing quantities, and 
we obtained simplified expressions for integer and half-integer 
values of n. These series expansions allow to perform arbitrary- 
precision calculations of the surface mass density and lensing 
properties of Einasto dark matter haloes. We also studied the 
asymptotic behaviour of the surface mass density, cumulative 
surface mass and deflection potential at small and large radii us- 
ing the series expansions. 

Furthermore, we compared the Sersic and Einasto surface 
mass densities using the equivalent values for the Sersic m and 
Einasto « indices for fixed values of T/e, Rg, por^ and r^, showing 



that both profiles have a similar behaviour. However, we noted 
that for the Einasto profile the external wings are more spread 
out and it seems to be less sensitive than the Sersic profile to the 
value of the surface mass density for a given Einasto index and 
radius in the inner region. Also, we find that the Einasto profile 
is more 'cuspy' than the Sersic profile: the former has higher 
values of the central surface mass density. These features are of 
key importance, because it is in this region that the lensing effect 
is more important, and these dissimilarities in the surface mass 
densities imply a difference in the lensing properties of both 
profiles. This result agrees with previous work of Cardone et al. 
(2005) and Dhar & Wifliams (2010). 

Our results are the first step in studying the properties of 
the Einasto profile using analytical means. The constant in- 
crease of computational power opens the possibility of using 
more realistic and sophisticated profiles like the Einasto pro- 
file in cosmological studies, where our results may apply. For 
example, they can be used in strong- and weak-lensing stud- 
ies of galaxies and clusters, where dark matter is believed to 
be the main mass component and the mass distribution can be 
assumed to be given by an Einasto profile. The better perfor- 
mance in cosmological A^-body simulations of the Einasto pro- 
file (Navarro et al. 2004; Merrittetal. 2006; Gao et al. 2008; 
Hayashi & White 2008; Stadel et al. 2009; Navarro et al. 2010) 
makes its inclusion in strong and weak lensing studies very 
promising. Recently, Chemin et al. (2011) analysed the rotation 
curves (RC) of spiral galaxies from THINGS (The HI Nearby 
Galaxy Survey, (de Blok et al. 2008)) and found that the Einasto 
profile provides a better match to the observed RC than the 
NEW profile (Navarro, Frenk, & White 1996, 1997) and the 
cored pseudo-isothermal profile. Also, Dhar & Williams (201 1) 
modelled the surface brightness profiles of a sample of elliptical 
galaxies of the Virgo cluster using a multi-component Einasto 
profile based on an analytical approximation, and obtained a 
good fit for shallow-cusp and steep-cusp galaxies with fit resid- 
ual errors lower in comparison to measurement errors in a wide 
dynamical radial range. Additionally, in their models of the most 
massive galaxies, the outer components are characterised by be- 
ing in the range 5 < « < 8, and a comparable range is ob- 
tained from A^-body simulations for the Einasto profile. Our ex- 
act analytical results for the spatial and lensing properties of the 
Einasto may be used to constrain the value of the Einasto in- 
dex and determine if the galaxy or cluster studied is dark matter 
dominated or not. More studies like Chemin et al. (2011) and 
Dhar & Williams (2011) could help to strengthen the position 
of the Einasto profile as a new standard model for dark matter 
haloes. Also, increasing the use of the Einasto profile in new 
cosmological studies could provide progress towards a solution 
to the cusp-core problem. 

This paper continues the effort initiated by Baes & Gentile 
(2011) and Baes & vanHese (2011) to advocate the use of the 
Fox H and Meijer G functions in theoretical astrophysics, in 
particular for studying the analytical properties of density mod- 
els like the Sersic model, where no additional analytical progress 
could be made until the Mellin integral transform formalism was 
applied. We hope that our work has again demonstrated the use- 
fulness of the Fox H and Meijer G functions as tools for analyt- 
ical work. 
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